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BEAM HARDENING POST-PROCESSING METHOD AND X-RAY CT 

APPARATUS 

BACKGROUND OF THE INVENTION 
5 The present invention relates to a beam hardening (sometimes abbreviated 

as 'BH' hereinbelow) post-processing method and an X-ray CT apparatus for 
correcting intensity of X-rays passing through a subject based on phantom data. 

An X-ray source employed in an X-ray CT apparatus outputs X-rays having 
a certain range of energy. On the other hand, the linear absorption coefficient 
10 for X-rays passing through a subject depends upon the X-ray energy, and 
exhibits a beam hardening effect by which the average energy shifts to a higher 
level as the transmission length through the subject becomes larger. Thus, the 
X-ray transmission intensity, i.e., projection information value, and the 
transmission length do not maintain a proportional relationship with each other, 
15 i.e., are non-linearly related to each other. 

The BH effect results in a cupping effect by which intensity decreases in the 
central portion of a reconstructed image, and is corrected by, for example, 
determining a corrective function for each channel of an X-ray detector for 
correcting the projection information values to produce a reconstructed image 
20 with homogeneous intensity (see Patent Document 1, for example). 

. Moreover, to achieve correction with higher accuracy, a plurality of 
cylindrical phantoms of different diameters, which diameters are sufficient to 
generally cover the entire FOV (i.e., imaged region), disposed at an imaging 
center are imaged, and projection information of the phantoms are used to 
25 improve the accuracy of the aforementioned correction. 

[Patent Document 1] 
Japanese* Patent Application Laid Open No. H5-130987 (pages 2 — 3, 
Figures 1 — 2). 

According to the conventional technique, however, a plurality of large and 
30 difficult-to-handle phantoms must be repeatedly imaged when acquiring 



projection information, which takes a long time; and at the same time, it is not 
possible to conduct correction with high accuracy in which a non-linear effect is 
factored into the projection information values. Specifically, since projection 
information value correction with high accuracy requires a large number of 
projection information values having different magnitudes for every channel, 
phantoms having a variety of diameters disposed at the imaging center must be 
imaged. 

Especially in acquiring calibration information for the X-ray CT apparatus, 
it is necessary to use two or three phantoms having diameters of 20 — 50 cm and 
conduct imaging that takes over 100 minutes only for improving the accuracy of 
the correction, which imposes heavy work on the operator. 

In consideration of such problems, it is important to find a way to 
implement a beam hardening post-processing method and an X-ray CT 
apparatus that can improve the accuracy of channel-by-channel correction on the 
BH effect easily and yet taking the non-linear effect into account. 

SUMMARY OF THE INVENTION 

Therefore, an object of the present invention is to provide a beam 
hardening post-processing method and an X-ray CT apparatus that can improve 
the accuracy of channel-by-channel correction on the BH effect easily and yet 
taking the non-linear effect into account. 

To solve the aforementioned problems and attain said object, a beam 
hardening post-processing method, in accordance with the invention of a first 
aspect, is characterized in comprising: obtaining one sinogram by acquiring first 
projection information of a phantom disposed at a position offset from an 
imaging center of an imaged region by imaging said phantom in a plurality of 
views from multiple directions; generating second projection information by 
performing beam hardening effect correction on said first projection information; 
generating third projection information by performing first function fitting on 
said second projection information; performing second function fitting on values. 
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of said third projection information using values of said second projection 
information in all said views as an independent variable for each channel 
constituting said second projection information; and correcting projection 
information of a subject disposed in said imaged region using a corrective 
5 function determined by said second function fitting. 

According to the invention of the first aspect, one sinogram is obtained by 
acquiring first projection information of a phantom disposed at a position offset 
from an imaging center of an imaged region by imaging the phantom in a 
plurality of views from multiple directions, second projection information is 

10 generated by performing beam hardening effect correction on the first projection 
information, third projection information is generated by performing first 
function fitting on the second projection information, second function fitting is 
further performed on values of the third projection information using values of 
the second projection information in all the views as an independent variable for 

15 each channel constituting the second projection information, and projection 
information of a subject disposed in the imaged region is corrected using a 
corrective function determined by the second function fitting; and therefore, 
when the corrective function is determined by the function fitting, the fitting is 
applied to the corrective function using a wide range of values of the second 

20 projection information because the values of the second projection information 
are different from view to view, and thus, the accuracy of the corrective function 
is improved. 

A beam hardening post-processing method, in accordance with the 
invention of a second aspect, is characterized in comprising: obtaining a plurality 

25 of sinograms by acquiring first projection information of a plurality of phantoms 
having different sizes disposed at a position offset from an imaging center of an 
imaged region by imaging each phantom in a plurality of views from multiple 
directions; generating second projection information by performing beam 
hardening effect correction on said first projection information; generating third 

30 projection information by performing first function fitting on said second 
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projection information; performing second function fitting on values of said third 
projection information using values of said second projection information in all 
said views and in all said sinograms as an independent variable for each channel 
constituting said second projection informaticm; and correcting projection 
5 information of a subject disposed in said imaged region using a corrective 
function determined by said second function fitting. 

According to the invention of the second aspect, a plurality of sinograms 
are obtained by acquiring first projection information of a plurality of phantoms 
having different sizes disposed at a position offset from an imaging center of an 

10 imaged region by imaging each phantom in a plurality of views from multiple 
directions, second projection information is generated by performing beam 
hardening effect correction on the first projection information, third projection 
information is generated by performing first function fitting on the second 
projection information, second function fitting is further performed on values of 

15 the third projection information using values of the second projection 
information in all the views and in all the sinograms as an independent variable 
for each channel constituting the second projection information, and projection 
information of a subject disposed in the imaged region is corrected using a 
corrective function determined by the second function fitting; and therefore, 

20 when the corrective function is determined by the second function fitting, the 
fitting is applied to the corrective function using a still wider range of values of 
the second projection information because the values of the second projection 
information are significantly different from sinogram to sinogram, and thus, the 
accuracy of the corrective function is further improved to improve image quality, 

25 or alternatively, the corrective function can be easily determined by obtaining the 
corrective function with high accuracy using a smaller number of phantoms. 

An X-ray CT apparatus, in accordance with the invention of a third aspect, 
is for acquiring projection information of an X-ray beam passing through an 
imaged region using an X-ray detector comprising multiple channels in a 

30 plurality of views from multiple directions, and performing beam hardening 
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effect correction on said projection information, and is characterized in 
comprising: obtaining means for obtaining one sinogram by acquiring first 
projection information of a phantom disposed at a position offset from an 
imaging center of said imaged region by imaging said phantom in all said views; 
5 beam hardening correcting means for generating second projection information 
by performing said beam hardening effect correction on said first projection 
information; first fitting means for generating third projection information by 
performing first function fitting on said second projection information; second 
fitting means for determining a corrective function by performing second 

10 function fitting on values of said third projection information using values of said 
second projection information in all said views as an independent variable for 
each channel constituting said second projection information; and correcting 
means for correcting projection information of a subject disposed in said imaged 
region using said corrective function. 

15 According to the invention of a third aspect, obtaining means obtains one 

sinogram by acquiring first projection information of a phantom disposed at a 
position offset from an imaging center of an imaged region by imaging the 
phantom in all views, beam hardening correcting means generates second 
projection information by performing beam hardening effect correction on the 

20 first projection information, first fitting means generates third projection, 
information by perforating first function fitting on the second projection 
information, second fitting means determines a corrective function by 
performing second function fitting on values of the third projection information 
using values of the second projection information in all the views as an 

25 independent variable for each channel constituting the second projection 
information, and correcting means corrects projection information of a subject 
disposed in the imaged region using the corrective function; and therefore, when 
the corrective function is determined by the second function fitting, the fitting is 
applied to the corrective function using a wide range of values of the second 

30 projection information because the values of the second projection information 
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are different from view to view, and thus, the accuracy of the corrective function 
is improved. 

The X-ray CT apparatus, in accordance with the invention of a fourth 
aspect, is characterized in that: said phantom has a circular cross-sectional shape. 
5 According to the invention of the fourth aspect, since the phantom has a 

circular cross-sectional shape, it is possible to change the value of the second 
projection information in a continuous manner. 

The X-ray CT apparatus, in accordance with the invention of a fifth aspect, 
is characterized in that: said cross-sectional shape has a diameter smaller than 
10 that of said imaged region. 

According to the invention of the fifth aspect, since the cross-sectional 
shape has a diameter smaller than that of the imaged region, work of the 
operator is mitigated by using a smaller phantom. 

The X-ray CT apparatus, in accordance with the invention of a sixth aspect, 
15 is characterized in that: said first function fitting is performed from channel to 
channel of said second projection information. 

According to the invention of the sixth aspect, since the first function fitting 
is performed from channel to channel of the second projection information, a 
corrective function can be determined in a channel direction. 
20 The X-ray CT apparatus, in accordance with the invention of a seventh 

aspect, is characterized in that: said first function fitting is performed from view 
to view of one sinogram containing a plurality of series of said second projection 
information. 

According to the invention of the seventh aspect, since the first function 
25 fitting is performed from view to view of one sinogram containing a plurality of 
series of the second projection information, smoothing is effected also in a view 
direction, and variation in the values of the second and third projection can be 
reduced when deteriruning the corrective function. 

The X-ray CT apparatus, in accordance with the invention of an eighth 
30 aspect, is characterized in that: said first function fitting comprises averaging 



means for averaging the values of said second projection information. 

According to the invention of the eighth aspect since the first function 
fitting averages the values of the second projection information by averaging 
means, smoothing can be easily achieved. 

i 

The X-ray CT apparatus, in accordance with the invention of a ninth aspect, 
is characterized in that: said first fitting means fits a higher-order function to the 
values of said second projection information. 

According to the invention of the ninth aspect, since the first function 
fitting fits a higher-order function to the values of the second projection 
information, when smoothing is performed in the view direction, for example, 
the smoothing can be achieved by the higher-order function fitting. 

The X-ray CT apparatus, in accordance with the invention of a tenth aspect, 
is characterized in that: said obtaining means obtains a plurality of said 
sinograms using a plurality of said phantoms having different diameters . 

According to the invention of the tenth aspect, since the obtaining means 
obtains a plurality of sinograms using a plurality of phantoms having different 
diameters, when the corrective function is determined by the function fitting, the 
fitting is applied to the corrective function using a still wider range of values of 
the second projection information because the values of the second projection 
information are significantly different from sinogram to sinogram, and thus/the 
accuracy of the corrective function is further improved to improve image quality, 
or alternatively, the corrective function can be easily determined by obtaining the 
corrective function with high accuracy using a smaller number of phantoms. 

The X-ray CT apparatus, in accordance with the invention of an eleventh 
aspect, is characterized in that: said beam hardening correcting means generates 
said second projection information by performing said beam hardening effect 
correction on said first projection information for each of said plurality of 
sinograms. 

According to the invention of the eleventh aspect, since the beam 
hardening correcting means generates the second projection information by 
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performing the beam hardening effect correction on the first projection 
information for each of the plurality of sinograms, the third projection 
. information can be generated for each sinogram. 

The X-ray CT apparatus, in accordance with the invention of a twelfth 
5 aspect, is characterized in that: said first function fitting generates said third 
projection information by performing said first function fitting on said second 
projection information for each of said plurality of sinograms. 

According to the invention of the twelfth aspect, since the first function 
fitting generates the third projection information by performing the first function 
10 fitting on the second projection information for each of the plurality of sinograms, 
a corrective f unction can be determined for each sinogram. 

The X-ray CT apparatus, in accordance with the invention of a thirteenth 
aspect, is characterized in that: said second fitting means determines said 
corrective function of a first order for each said sinogram. 
15 According to the invention of the thirteenth aspect, since the second fitting 

means determines a first-order corrective function for each sinogram, calculation 
efficiency is improved by determining approximate corrective functions. 

The X-ray CT apparatus, in accordance with the invention of a fourteenth 
aspect, is characterized in that: said second fitting means determines a 
20 higher-order corrective function by performing higher-order function fitting on 
said corrective functions of a first order. 

According to the invention of the fourteenth aspect, since the second fitting 
means determines a higher-order corrective function by performing higher-order 
function fitting on the corrective functions of a first order, a higher-order 
25 corrective function containing a non-linear effect can be easily determined. 

The X-ray CT apparatus, in accordance with the invention of a fifteenth 
aspect, is characterized in that: said second fitting means determines the 
corrective function by performing function fitting on the values of said third 
projection information for all said sinograms. 
30 According to the invention of the fifteenth aspect, since the second fitting 
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means determines the corrective function by performing function fitting on the 
values of the third projection information for all the sinograms, a corrective 
function with high accuracy can be determined from a large amount of the 
second and third projection information. 
5 The X-ray CT apparatus, in accordance with the invention of a sixteenth 

aspect, is characterized in that: said corrective function contains higher-order 
terms. 

According to the invention of the sixteenth aspect, since the corrective 
function contains higher-order terms, a non-linear effect can be included to 

10 achieve fitting with high accuracy 

According to the present invention, obtaining means obtains one sinogram 
or a plurality of sinograms by acquiring first projection information of a phantom 
or a plurality of phantoms of different diameters disposed at a position offset 
from an imaging center of an imaged region by imaging (each phantom in all 

15 views, beam hardening correcting means generates second projection 
information by performing beam hardening effect correction on the first 
projection information, first fitting means generates third projection information 
by performing first function fitting on the second projection information, second 
fitting means determines a corrective function by performing second function 

20 fitting on values of the third projection information using values of the second 
projection information in all the views as an independent variable for each 
channel constituting the second projection information, and correcting means 
corrects projection information of a subject disposed in the imaged region using 
the corrective function; and therefore, when the corrective function is determined 

25 by the function fitting, the fitting is applied to the corrective function using a 
wide range of values of the second projection information because the values of 
the second projection information are different from view to view or from 
sinogram to sinogram, and thus, the accuracy of the corrective function is 
improved to improve image quality, or alternatively, the corrective function can 

30 be easily determined by obtaining the corrective function with high accuracy 
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using a smaller number of phantoms. 

Further objects and advantages of the present invention will be apparent 
from the following description of the preferred embodiments of the invention as 
illustrated in the accompanying drawings. 

5 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is a block diagram showing the general configuration of an X-ray 
CT apparatus. 

Figure 2 is a functional block diagram showing a data processing apparatus 
10 of Embodiment 1 . 

Figure 3 shows a positional relationship between a phantom and a rotating 
section of Embodiment 1. 

Figure 4 shows a sinogram and projection information values of the 
phantom of Embodiment 1. 
15 Figure 5 is a flow chart showing the operation of the data processing 

apparatus of Embodiment 1 . 

Figure 6 is a block diagram showing files in a storage device of 
Embodiment 1. 

Figure 7 shows processing on projection information values in a channel 
20 direction in accordance with one embodiment. 

Figure 8 shows processing on the projection information values in a view 
direction in accordance with the embodiment. 

Figure 9 shows correction factors for the projection information values in 
accordance with the embodiment. 
25 Figure 10 is a diagram explaining determination of a second fitting function 

for the projection information values in accordance with the embodiment. 

DETAILED DESCRIPTION OF THE INVENTION 

Preferred embodiments of a beam hardening post-processing method and 
30 an X-ray CT apparatus in accordance with the present invention will now be 
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described with reference to the accompanying drawings. 

The general configuration of the X-ray CT apparatus of the present 
embodiment will be first described. Figure 1 shows a block diagram of an X-ray 
CT apparatus. As shown in Figure 1, the apparatus, comprises a scan gantry 2 
5 and an operating console 6. 

The scan gantry 2 has an X-ray tube 20. X-rays (not shown) emitted from 
the X-ray tube 20 are formed into a fan-shaped X-ray beam, i.e., fan-beam X-rays, 
for example, by a collimator 22, and projected toward an X-ray detector 24. 

The X-ray detector 24 has a plurality of X-ray detector elements arranged in 
10 line as an array in the extent direction of the fan-beam X-rays. The X-ray 
detector 24 is a multi-channel detector in which the plurality of X-ray detector 
elements are arranged in line as an array. 

The X-ray detector 24 as a whole forms an X-ray impingement surface 
curved as a cylindrical concavity. The X-ray detector 24 is formed of a 
15 combination of scintillators and photodiodes, for example. It should be noted 
that the X-ray detector 24 is not limited thereto but may comprise semiconductor 
X-ray detector elements using cadmium telluride (CdTe) or the like, or ionization 
chamber X-ray detector elements using xenon (Xe) gas, for example. The X-ray 
tube 20, collimator 22 and X-ray detector 24 together constitute an X-ray 
20 emitting/detecting apparatus. 

The X-ray detector 24 is connected with a data collecting section 26. The 
data collecting section 26 collects data detected by the individual detector 
elements in the X-ray detector 24. The emission of the X-rays from the X-ray 
tube 20 is controlled by an X-ray controller 28. The interconnection between the 
25 X-ray tube 20 and X-ray controller 28 and that between the collimator 22 and a 
collimator controller 30 are omitted in the drawing. The collimator 22 is 
controlled by the collimator controller 30. 

The above-described components from the X-ray tube 20 through the 
collimator controller 30 are mounted on a rotating section 34 of the scan gantry 2. 
30 In this configuration, a subject or a phantom is laid on a cradle within a bore 29 in 
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the center of the rotating section 34. The rotating section 34 is controlled by a 
rotation controller 36 to rotate, while X-rays are emitted from the X-ray tube 20 
and X-rays passing through the subject or phantom are detected at the X-ray 
detector 24 as projection information for each view. The interconnection 
5 between the rotating section 34 and rotation controller 36 is omitted in the 
drawing. 

The operating console 6 has a data processing apparatus 60. The data 
processing apparatus 60 comprises, for example, a computer. The data 
processing apparatus 60 is connected with a control interface 62. The control 

10 interface 62 is connected with the scan gantry 2. The data processing apparatus 
60 controls the scan gantry 2 via the control interface 62. 

The data collecting section 26, X-ray controller 28, collimator controller 30 
and rotation controller 36 in the scan gantry 2 are controlled via the control 
interface 62. The individual connections between these sections and the control 

15 interface 62 are omitted in the drawing. 

The data processing apparatus 60 is also connected with a data collection 
buffer 64. The data collection buffer 64 is connected with the data collecting 
section 26 in the scan gantry 2. Data collected at the data collecting section 26 
are input to the data processing apparatus 60 via the data collection buffer 64. 

20 The data processing apparatus 60 perforins image reconstruction using 

transmitted X-ray signals, or projection information, collected via the data 
collection buffer 64. The data processing apparatus 60 is also connected with a 
storage device 66. The storage device 66 stores projection information collected 
in the data collection buffer 64, reconstructed tomographic image information, 

25 programs for implementing the functions of the present apparatus, and so forth. 

The data processing apparatus 60 is further connected with a display 
device 68 and an operating device 70. The display device 68 displays the 
tomographic image information and other information output from the data 
processing apparatus 60. The operating device 70 is operated by an operator, 

30 and supplies several kinds of instructions and information to the data processing 
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I. 

apparatus 60. The operator interactively operates the present apparatus using 
the display device 68 and operating device 70. The scan gantry 2, an imaging 
table 4 and the operating console 6 together constitute an acquiring apparatus for 
imaging the subject or phantom to acquire a tomographic image. 
5 Figure 2 shows a functional block diagram of only a portion associated 

with the beam hardening post-processing method, which is one embodiment of 
the present invention, in the data processing apparatus 60. The data processing 
apparatus 60 comprises the functions of BH correcting means 201, first fitting 
means 202, second fitting means 204 and image reconstructing means 206 

10 performed on the projection information in the storage device 66. 

The BH correcting means 201 performs BH correction on the projection 
information in the storage device 66. Representing a value of the projection 
information acquired at each channel in the X-ray detector 24 as Ih, and BH 
corrected data as Ic, the BH correction is conducted according to the following 

15 equation: 

Ic = Bo • Ih + Bi • + B 2 - Ih 3 + B 3 • Ih*, (1) 

where B0 — B3 are correction factors. The correction factors are established for 
each channel by a method described in Patent Document 1, for example, and are 
saved in the storage device 66 in a correction factor table. 

20 The first fitting means 202 performs channel-to-channel or view-to-view 

smoothing on the projection information in the storage device 66. The first 
fitting means 202 conducts the fitting by averaging means for performing 
channel-to-channel or view-to-view average calculation, or by fitting a 
higher-order function to the projection information values in a channel or view 

25 direction. A function determined by the fitting offers a similar effect to that by 
smoothing because high frequency components beyond the order of the function 
are removed. 

The second fitting means 204 performs first-order or higher-order function 
fitting on projection information values acquired by one channel in the X-ray 
30 detector 24 and subjected to the first function fitting by the first fitting means 202. 
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Thus, a corrective function similar to Equation (1) used by the BH correcting 
means 201 can be obtained. 

The image reconstructing means 206 reconstructs a tomographic image of 
the subject or phantom using a sinogram that is projection information 
comprised of a plurality of views in the storage device 66. The image 
reconstruction is conducted by using a filtered backprojection method, for 
example, and the reconstructed image is displayed on the display device 68. 

Next, before explaining the operation of the X-ray CT apparatus in 
accordance with one embodiment of the present invention, data collection, 
projection information and a sinogram for a phantom disposed at a position 
offset from the imaging center of the bore 29 will be described. 

Figure 3 shows a phantom 310 disposed in the bore 29 of the scan gantry 2. 
The phantom 310 has a circular cross section, and its center lies at a position 
different from the imaging center of the bore 29. An X-ray fan beam generated 
from the X-ray tube 20 and passing through a bow-tie filter passes through the 
phantom 310, and is detected by the X-ray detector 24. 

The X-ray detector 24 in which a plurality of X-ray detector elements are 
arranged in line in a array in the extent direction of the X-ray fan beam detects 
projection information of the phantom 310 at individual channels in the array. 
The X-ray detector 20, collimator 22 and X-ray detector 24 are opposingly 
disposed with respect to the bore 29, and conducts projection information 
acquisition while rotating around the bore 29 along with the rotating section 34 
without changing their relative positions. The projection information is 
acquired for each view index that corresponds to a rotation angle, and one 
sinogram is generated. 

Figure 4(A) shows an exemplary sinogram when the phantom 310 is used. 
The sinogram is comprised of a projection information part lying near the center 
of the channels and an air data part lying near the peripheries of the channels. 
Since the phantom 310 is disposed offset from the imaging center, the position of 
a width of channels of the projection information part varies in position with the 
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rotation of the rotating section 34, i.e., with the varying view index, and 
meanders in the view index direction as shown in Figure 4(A). For the same 
reason, the width of channels of the projection information part moves with the 
varying view index. 

Figure 4(B) is a diagram in which projection information at a view index of 
in Figure 4(A) is represented with a horizontal axis of the channel index and a 
vertical axis of the projection information value. Since the projection 
information value is proportional to the transmission length of the X-ray beam 
passing through the phantom 310, X-rays passing through near the center of the 
phantom 310 exhibit a higher projection information value because of a larger 
transmission length, and those passing through near the peripheries of the 
phantom 310 exhibit a lower projection information value because of a smaller 
transmission length, resulting in a semicircular projection image as shown in 
Figure 4(B). 

As an example, consider a projection information value at a view index of j 
and a channel index of i. An X-ray beam indicated by a broken line in Figure 3 
impinges upon the channel index i in the X-ray detector 24 when the vieiw index 
is j. Define the length over which the X-ray beam passes through the phantom 
310 at that time as Z. The relationship of the length I to a projection information 
value h at the channel i in Figure 4(B) is: 

loch. . 
Moreover, since the phantom 310 lies offset from the imaged region in Figure 3, 
the transmission length I of the channel i varies from view to view. Therefore, 
the projection information value h of the channel i shown in Figure 4(B) varies 
from view to view. 

Figure 4(C) is a diagram in which the projection information value of a 
channel index of i in Figure 4(A) is represented with a horizontal axis of the view 
index and a vertical axis of the projection information value. Since the 
projection information value is proportional to the transmission length of the 
X-ray beam passing through the phantom 310, which length varies from view 
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index to view index, a periodic function as shown in Figure 4(C) results. 

Next, the operation for determining the correction factors in the BH 
post-processing will be described. Figure 5 is a flow chart showing the 
operation for determining the correction factors in accordance with an 
embodiment of the present embodiment. Moreover, Figure 6 shows files of 
intermediate projection information created in this operation. A phantom is 
first disposed at a position offset from the imaging center in the bore 29. The 
phantom is made from a material such as polypropylene, and has a cylindrical 
shape with a diameter of 35 cm, for example. The phantom is used to conduct a 
phantom scan (Step S501). In Figure 6, first projection information 601 acquired 
by the scan is shown. A sinogram comprised of the first projection information 
is subjected to pre-processing such as noise removal and sensitivity correction 
(StepS502). 

, Then, BH correction is applied to a projection information value Ih using 
Equation (1) to obtain a corrected projection information value Ic (Step S503). 
Second projection information shown in Figure 6 is thus generated. In this file, 
the BH effect is largely removed, but some BH effect due to channel-by-channel 
variation in the X-ray detector 24 is still left. Exemplary second projection 
information is schematically shown in Figure 7(A). Although the second 
projection information has a generally semicircular shape that represents 
projection information of the circular phantom, the projection information value 
Ic varies in a pulse-like manner because of differences in X-ray sensitivity among 
some channels, for example. The variation must be corrected channel by 
channel because it is a phenomenon specific to individual channels. Moreover, 
Figure 8(A) schematically shows exemplary projection information values in the 
second projection information of one channel in the view direction. The 
projection information value Ic varies in a pulse-like manner among some view 
indices. 

Referring again to Figure 5, the second projection information 602 is then 
used to perform smoothing in the channel direction by the first fitting means 202 
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(Step S504). This step generates third-prime projection information 603 shown 
in Figure 6. In the projection information, the projection information value Ic 
resulting from the channel-to-channel variation is smoothed and removed. 
Exemplary third-prime projection information is schematically shown in Figure 
7(B). Only, the semicircular shape that represents projection information of the 
circular phantom is obtained as the projection information. 

Referring again to Figure 5, the third-prime projection information 603 is 
then used to perform smoothing in the view direction by the first fitting means 
202 (Step S505). This step generates third projection information 604 shown in 
Figure 6. In the projection information, the projection information values Ic 
resulting from the view-to-view variation in one channel is smoothed. 
Exemplary third projection information is schematically shown in Figure 8(B). 
The periodic projection information values of one channel in the view direction 
are smoothed. 

Referring again to Figure 5, a first-order corrective function is determined 
from the second and third projection information using the second fitting means 
204 (Step S506). Representing a projection information value of the second 
projection information at a channel index of i as S(j), a corresponding projection 
information value of the third projection information as F(/), and plotting the 
projection information values for all view indices with a horizontal axis of S(/) 
and a vertical axis of F(/), the projection information values are lined up on a 
straight line generally passing through the origin, as shown in Figure 9. The 
straight line is defined as a corrective function for the channel i. Moreover, the 
corrective function is saved as corrective function information 605 in the storage 
device 66. Representing the gradient of the corrective function as Ki, the 
following relationship holds: 
F(/)/S(/>Ki ; 

Multiplying a BH-corrected projection information value Ic of the channel i 
acquired from a subject by the correction factor Ki gives: 
Ip = Ic - Ki, 
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and a smoothed and corrected projection information value Ip of the subject is 
thus obtained. 

As shown in Figure 9, magnitude of the projection information value S(j) 
and the acquired value range depend upon the diameter and the position in the 
5 bore 29 of the phantom 310 because the magnitude of the projection information 
value is proportional to the transmission length I shown in Figure 3. 

Referring again to Figure 5, whether or not to improve the accuracy of the 
corrective function is determined (Step S507). If the accuracy of the corrective 
function is to be improved (Yes in Step S507), a phantom of a diameter different 
10 from that of the phantom 310 is disposed at a position different from the imaging 
center in the bore 29, and the flow goes back to Step S501. Then, a new 
corrective function is obtained. 

Figure 9(B) shows exemplary corrective functions obtained by using two 
phantoms of different diameters. A transmission length of an X-ray beam 4 
15 depends upon the diameter and position of the phantom, and so does a 
projection information value S(/). Therefore, when the diameters of the 
phantoms are represented as A and B, and if A < B, the projection information 
values of the phantom A lie in a region A and those of the phantom B lie in a 
region B, as generally depicted in Figure 9(B). The corrective functions are 
20 determined from the projection information values in these regions. 

Thereafter, if sufficient data in view of the accuracy of the corrective 
function has been acquired (No in Step S507), a higher-order function fitting is 
applied to the correction factors in each of the acquired regions (Step S508). 
Figure lO shows a case in which the corrective functions for the phantoms A and 
25 B shown in Figure 9(B) are employed. A third-order fitting function given 
below is fit to values of the corrective function A in the region A and of the 
corrective function B in the region B: 

If = K0 - S(j) + Kl - S(/)2 + K2 • S(/) 3 , (2) 

to determine the correction factors K0, Kl and K2. At that time, the correction 
30 factors in the region A having smaller projection information values may be 
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regarded as having higher accuracy than those in the region B having larger 
projection information values, and therefore, it is possible to apply weighting for 
each region and determine the correction factors in Equation (2) so that the fitting 
is effected more accurately in the region A. 
5 Referring again to Figure 5, higher-order corrective function information 

606 comprised of the correction factor values K0, Kl and K2 is then saved in the 
storage device 66 (Step S509), and the process is terminated. 

When conducting imaging of a subject, the correction factors K0, Kl and K2 
for each channel are applied to BH-corrected projection information values Ic of 
10 the subject to determine corrected projection information values If from Equation 
(2). The projection information values If are then subjected to image 
reconstruction by the image reconstructing means 206 to obtain tomographic 
image information. 

As explained in the preceding description, in the present embodiment , 

15 since phantoms of different diameters are disposed at a position offset from the 
imaging center and projection information having the transmission length of the 
X-ray beam varying from view to view, and hence the projection information 
value varying from view to view, is acquired for each channel, 
channel-by-channel correction on the projection information values conducted 

20 after BH correction can be achieved by approximating a corrective function using 
a higher-order function while taking non-linear components into account, and 
moreover, a corrective function with high accuracy can be determined using a 
smaller amount of phantom data, therefore mitigating physical work and 
reducing the time for the work for calibration by the operator. 

25 Moreover, although the fitting is performed using a third-order function 

using Equation (2) in this embodiment, the fitting may be performed using 
second-order or fourth- or higher-order function. 

Furthermore, although a first-order corrective function is determined for 
each sinogram at Step S506 in the embodiment, the second projection information 

30 602 and third projection information 604 for each sinogram are fitted with a 



-20- 



higher-order function to determine a higher-order corrective function without 
determining the first-order corrective function. 

Many widely different embodiments of the invention may be configured 
without departing from the spirit and the scope of the present invention. It 
should be understood that the present invention is not limited to the specific 
embodiments described in the specification, except as defined in the appended 
claims. 



